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Analytic dynamics of one-dimensional particle system 

with strong interaction 

V. A. Malyshev 

Abstract 



We study here the small time dynamics of N electrons on the circle with Coul;omb 
repulsive interaction and study the series for the velocities (initially zero) . The main result 
is the estimates of the convergence radius from below. We explain how this result is related 
Q-ij to the problem of very fast establishing of direct electric current. 



1 Introduction 



Model and results We consider a system of N point particles i — 1, 2, ..., JV on the interval 
>• ! [0, L) & R, with periodic boundary conditions, that is on the circle Sl of length L. Initially, 

^ \ they are situated at the points 

00 . 

= xi(0) < ... <x N (0) < L 

en 

O ■ The trajectories Xi(t) are defined by the following system of N equations 



d 2 Xi dU 

-w = -d X - + F{x * ) (1) 

where the interaction between the particles is 

where summations is over all pairs of nearest neighbors on the circle. The Coulomb potential 

V{x) = V(—x) — -,r — \x\ is assumed, and we denote f(r) = -ijp- = r~ 2 the interaction 

force. Note that the potential is repulsive, infinite at zero and thus, the particles cannot change 
their order . F(x) is the external force. 

In the papers [TT| [T2] we considered fixed points of such systems. Dynamics is more com- 
plicated to understand. It is standard that the solution of the system ([T|) exists on all time 
interval and is unique (for any initial conditions), however to get more detailed information 
about particle trajectories (for sufficiently large N) is sufficiently difficult and demands elabo- 
ration of special methods. If moreover F is analytic, then it is well-known [lj, that the solution 
can be presented as the convergent power series in t in some neighborhood of t = 0. 

Here we consider natural initial conditions, that is for all i 

A = A,(0) = x l+1 (0) - Xi (0) = ^,Vi(0) = (2) 



and moreover, it is convenient to put £i(0) = 0. Note that this configuration is a fixed point 
for zero external force. We are looking for the solution of the form 

oo 

3=1 

and get bounds for the convergence radius, dependent on N, of these series under the initial 
conditions (|21). 

Theorem 1 Let F be analytic on the circle Sl- Then 

1. for j = 1,2,..., there exist numbers bj < oo, not depending on N and such, that for all 
i,j and any N 

\cij\ < bjN 2 

2. if moreover, for some constant Cp > and all x and k 

\F {k \x)\ <C k F + \ 
then there exists constant < x < °°; n °t depending on N and such that for all i,j 

5 ■ 3 

\cij\ < x^Ns 3 '? 

It follows that the convergence radius R = R(N) of the series (J3J) has a bound from below 
R > x _1 iV - 6. From the proof of the first assertion of the theorem (see section 2.2) one can 
see that the bound from above could be of the order -j=, but this not yet rigorously proved. 
During the proof of the second assertion of the theorem we give the explicit bound for x- Also, 
explicit formulas for Cij,j = 1, 2, 3, 4 are presented. 

Why this model Mathematical problems of statistical physics are elaborated sufficiently 
deeply for equilibrium systems on the lattice. But on the continuous (euclidean) space only for 
particle systems with small inverse temperature or small density. There are many other cases 
where the problems are not even formulated on the mathematical level. One of such cases is 
the direct electric current. It is described by Ohm's law on the macro level. On the micro level, 
in any textbook on condensed matter physics, it is described as a flow of free (or almost free) 
electrons, any of which is accelerated by the external force and impeded by the external media 
(crystal lattice). Physics and mathematics introduced a lot of such one-particle models with 
constant accelerating force and various models (about 20, historically the first such model is 
the famous Drude's model of 1900) of external media, where the particles loose their kinetic 
energy. 

The central question is where the accelerating force comes from. The problem is that the 
power lines can have hundreds kilometers of length but the external force acts only on the 
length of several meters of the wire. In fact, there are even more problems with the direct 
electric current. For example, one should explain why the speed of the flow is permanent and 
sufficiently slow, but this regime is being established almost immediately. We discussed the 
first problem in (|10j). but there was no Coulomb interaction there. We shall come back to this 
problem in the next paper. Here we discuss the second problem. 

Consider first a trivial case with constant F > and with initial conditions (j2j). Then it is 
clear that for any % 

Ft 2 

v i (t) = Ft,x t {t)=x i {0) + — 



that is Xi(t) are analytic in t for any t. Moreover, to reach the speed of order 1 (for example 
Vi(t) = 1) the time of order 1 is necessary. We shall see however that the situation is completely 
different for non-constant F, Anticipating the events, it is important to note that technically 
the difference occurs because for constant F all discrete derivatives in the recurrent equations 
below are identically zero. 

If we could prove that the coefficients Qj in the expansion Vi = J2^=i c i,j& grow as N a i for 
some a > 0, then the speed of order 1 will be achieved much earlier - for time of the order 
t = N~ a , that is almost immediately. We cannot prove such result as it is but the estimates of 
the coefficients, obtained here, make such result quite plausible. 

The techniques of the paper is apparently new. There are results for small time dynamics 
of multi-particle systems (see, [2]-|9]) but they are not applicable in our case because of very 
strong interaction. 

2 Proof 

2.1 Equations for the coefficients 

Fix initial data Xi(0),Vi(0) as in (|2j) and consider the trajectories Xi(t) G Sl on the interval 
< t < t for some t = t (N) > 0. Putting 

Ai(t) = x l+1 (t) - Xi (t),A = Ai(0) = — 
we have the equations 

—^ = f{xi{t) - x^)) - f(x i+ i(t) - Xi (t)) + F{ Xi (t)) 
at 

or 

-^ = f(A+l [v i (t 1 )-v i - 1 {t 1 )]dh)-f{A+ I [v l+1 (t 1 )-v l (t 1 )})dt 1 )+ 

+F(xi(0) + I Vi(ti)dti) 
Jo 

Integral equations Equivalent system of integral equations is 

Vi {t) = [ [/(A+ f [viih) -«i_i(ti)]dti)- 



- /(A + / [v i+1 {h) - u«(ti)]d*i) + F{ Xi {0) + f Viitjdt^dt (4) 

Jo Jo 

and can be rewritten as follows 

Vi {t) = I ((A + Ri-tit))- 2 - (A + R t (t))- 2 + F{xi{0) + f v 4 (ti)dti))dt (5) 

Jo Jo 

where 

ft 



Ri-i(t)= / (viih) - u<_i(ti))dti 



For the sequel we need some notation for discrete derivatives. Let on the interval [0, N] C Z 
with periodic boundary conditions a function g(i) be given (that is a periodic function on Z 
with period N). Let us call 

(Vg)(i) = (V + g)(i) = g(i + 1) - g{{), (V"^)(z) = g(i) - g(i - 1) (6) 

its right and left derivative correspondingly. Note that they commute and the Leibnitz formula 
holds 

V+(gf)(i) = f(i + l)(V + </)(i) + <?«(V + /)« = (Sf)(V + g) + g(V + f) (7) 

where S is the shift operator 

(5/)(<) = /(» + l) 

Below discrete derivatives will act on the indices i. If the function f(i) does not depend on i, 
then its (discrete) differentiation gives zero. 

We come back to the main equations and rewrite them as follows 

Vi (t)= [ dt[(-\/-((A + R i (t))- 2 ) + F(x i (0)+ [ Vi(ti)dti)] 
Jo Jo 

The following representation of the integrand will be useful 

(A + Ri-^t))- 2 - (A + R t (t))- 2 + F(a*(0) + f v^dtx) = 



o 



= A- 2 (l + ^i)- 2 - A" 2 (l + ^)- 2 + F(a*(0) + f «,(«!)*!) = 

oo rt 

= F{ Xi {0)) + £ rf m A- 2 - m (i?ri - RT) + [i^MO) + / Wi(ti)dti) - F(x,(0))] 

771=1 •'O 

where 

d m = (-l) m (m + l) 

If F is analytic on Sl, then there exists e > sufficiently small and such that for any xq G Sl 
the following series 

F(x) = F(x ) + f] ^^(x - x ) k 
converges for all a; G [a?o — e, xo + e). Then finally 

««(*) = F( Xi (0))t + ft d m A- 2 - m [-V-RT)]dt + £ f F(fc) ^(0))(io^(^i)^i) fc dt (g) 

>/0 771=1 fc=1^0 fc - 

Recurrent equations Using ([3]) and 

Ri-l(t) = J2(°i,J ~ C i-hj)~— T ( 9 ) 

i=i ■> + l 

oo £j+l 

Fi - i?i_i = E(c i+ ij - 2qj + Ci-ij)— — (10) 

i=i -^ + - 1 

we see, substituting (jH|) to (ISlh that the right-hand side of flSJ) can also be presented as a 
convergent series with well-defined coefficients. 



We shall find Qj by equating the coefficients of P. For j = 1,2 the equations give immedi- 
ately 

c il = F(x i (0)),c i)2 = (11) 

as other summands in the right side of (jHJ) have larger order in t. For j > 3 the equations for 
the coefficients of P are 

J m=l fc=l rc! 

where, for the power series </>(£) = Y^k=o a kt k , we denote [</>(£)]j = %. For j > 2 the coefficients 
Cjj can be found recursively, moreover Qj depend only on q^ with k < j — 2. In fact, the 
right-hand part of the equation for Cjj does not contain c^ with A; > j — 1, as due to Qj, each 
of Cjfc appears together with t fe+1 . 
Then the main equations will be 



Cij 



- £ d m A- 2 -(-V-[(E(^ + ij - c„)-— ) m ],-i) + (13) 

J m=l i=l J i l 



m=l j=l 



FW(x,(0)) r ^ P+ 1 



+£^n^X>^- 



We have 



r(y Cjj i — )Vi= y _^___^_ (14 ) 

where J2j 1 +...+j rn =j- m -i is the sum over all ordered arrays ji, ...,jk, such that 

(ii + 1) + ... + (i fe + 1) = fe + ji + ... + j k = j-i (is) 

It follows 

k<Ji + ...+J k =J-l-k<j-2,k< l 3 -^] (16) 

Similarly 

i=i •/ + x h,-Jm Jl + - 1 - ■ /m + i 

and both (1T51) and ( !T6l) hold with & instead of m. 
That is why the equations can be written as 

c = Gc + c (o) (17) 

where c is the vector c = {c^}, the free term c^ = {cL } is simple 

cS 1 ) = F(x J (0)),4 0) = 0,j>2 (18) 

and non-linear operator G is defined by 

°a = c a = F(xi(Q)),c i2 = 



t¥i t¥i 



Cij = {Gc)ij = - Y, Yl Aj(m;ji,...,j m )+J2 Yl B ij (k',j 1 ,...,j k ) (19) 

m=l ii+...+j m =j'-m-l fe=l J!+...+j fe =j-fc-l 

for j > 3, where 

M^'Ju-M = -d m A- 2 - m v-(^f...^^) (20) 

J jl+1 Jm + 1 

Further on, i*i,jfe,g will denote any discrete derivative (Ilp=i V s ^)-F^(xj(0)), where s(p) = ±. 
For the estimates the choice of s(p) does not matter. Put F^k = -Fj^o- 

Explicit expression for c^, c^ is immediately obtained if in the equations (TT9I) we take into 
account only the terms with k — 1 and m = 1, as k.m < [^y-] < 1) 

c i3 = -^!A- 3 V-V + ^ + ^F«(x,)^ = kdiA-'F^ + Fi.o.oFi.i.o) 
o 2 o 2 o 

<* = -^A- 3 V-(V + f ) 2 + V)(* 4 4 = 

= -(-diA~ Fj A 2-Fi,o ; i + 2^,1,0^,0,0) 
From the formulas 

^,0,1 = ^+1,0,0 - ^,0,0 = f ' F^{x)dx, |F,, ,i| < C^A (22) 

f X ' _1_ 1 1" 'V I /\ 

^,0,2 = (#+2,0,0 - #+1,0,0) - (#+1,0,0 - #,0,0) = / ( / F i2 \y)dy)dx, |# A2 | < C 3 F A 2 

J X; J X 



it follows that 



|<te| < 3CKA- 1 + -), | Q4 | < -C| + -C\ 



It is easy to see that for most i the coefficients 0^3 are really of the order A l . 

2.2 Principal exponent 

From the recurrent formulas (120]) and (|21|) it follows that Cjj are finite and depend on i, j, AT. 
First of all, we shall study them as functions of N for fixed i,j. Otherwise speaking, we shall 
prove the first part of the theorem. Define the principal exponent 

I(£) = hmsup — - 

for a variable £ depending on N. Roughly speaking, it shows that the main order of the 
asymptotics of £ is N I( -&. 



We shall consider the algebra A of polynomials of countable number of (commuting) vari- 
ables Fi t k,q,i = 1) ■■■,N;k,q = 0,1,2,... with real coefficients, not depending on F. For any 
monomial M of this algebra denote 

Q(M) = - y £q 

over all q in this monomial. The natural mapping of A onto its subalgebra A , generated by 
all Fik = Fjfco, is defined by the subsequent substitutions 

"i,k,q "i-\-\,k,q— 1 "i,k,q— 1 

or with the following formula 

n 

(v+) n = (s - i) n = Y, c*(-i) k s n ~ k 

fc=0 

Lemma 1 For any monomial MgA 

I(M) < Q{M) 

It is sufficient to prove that 

I{Fi, q ) < Q(F ijq ) = -q 

This can be done as in (j22|) . using induction in q. Put 

gi, n = V n gi,A = — 
Then 

rXi+A ryi+A ry n -i+A 

g it n + i = V n+1 gi = g i+1 , n - g^ n = dy 1 {l dy 2 ....( dy n g {n) (y n ))) 

JXi J J/i Jy-n-i 

and if \g( n \x)\ < C™ +1 , we have 

h, n \ < c; +1 a" 

For any polynomial P = J2 a r M r with (different) monomials M r and coefficients a r , not 
depending on F, but possibly depending on N, define 

Q(P) = max(/(a r ) + Q{M r )) 

in agreement with the previous definition. Then for any polynomial P 

I(P) < max(J(or) + I(M r )) < max(I(a r ) + Q(M r )) 

Note that for any two polyonomials 

Q(PiP2)<Q(Pi) + Q(P2) 

Also 

Q(V+P) < Q(P) - 1, Q(V~V + F) < Q(P) - 2 (23) 

Denote degP the degree of the polynomial P = J2 a rM r , that is the maximal degree of its 
monomials. 



Lemma 2 For j > 2 Cij is a polynomial in the algebra Ao and has degree not greater than 

We already saw this for j = 3,4. Note that 

deg(V ± P) = degP 

Now one can use induction: in the formula ( 120]) the degree is j — 2, and in (}2~Tj) the degree will 
be j - 1. 

Note that the recurrent formulas define c^ for all functions F(x), not necessary analytic. 
That is why the following assertion makes sense. 

Lemma 3 Let F be infinitely differentiable. Then for all i,j 

j-1 

I(Cij) < QiCij) < — ^— 

As 

Q( Cij ) = 0, j = 1,2,Q(q i3 ) = l,Q(ci >4 ) = 

then the assertion holds for j = 1,2,3,4. We shall prove the lemma by induction in j. Assume 

that 

j-1 

Q{cij) < ~y~ 

forall j = 1,2,..., J -2. 

Then for given m,ji, ■■■,j m , accordingly to ( 1151) and (1161) . we have 



Q(Aij(m; j 1} ...,j m )) < 2 + m - 1 + Q{c ih ) + ... + Q{c ijm ) - m < 1 + -(ji + ... + j m ) - — = 

1 . m 

= 1 + - (J-m-1) 

2 V ; 2 

as, according to ( 12"3"|) . (—1) and (— m) are appended because of the discrete differentiation of 

the corresponding monomials. The last expression attains its maximum when m = 1. It follows 

Q(Aij(m;j u ...,j m )) < ~^— 

Similarly, for Bij(k : j\, ...,jk) we have the following inequalities 

Q(B u (k; 3l , ...,j k )) < i(J - 1 - k) - \ < J -=^- 

We get thus Q(aj) < ^-, and /(qj) < Q(qj) < ^. 

2.3 Convergence radius 

Here we shall prove the second assertion of the theorem 1. In the proof it is convenient to write 
TV instead of j; an d assume that Cp > 1. 

We shall use the following majorization principle for infinite system of recurrent equations 
and inequalities: for example, if two systems of equations are given 

cg ) =P®(cg ) ) ..-,cg_ a ),«=l,2 

where P^ are the polynomials with coefficients p^\ where p^ > 0, \p$\ < p^ for all a, and 
also \c\j | < c\j for j = 1,2,3,4, then |cL | < cL for all j. One of such system (we use it in the 
proof) corresponds to one-particle problem (that is with N = 1) with specially chosen external 
force, will be now introduced. Other auxiliary system /?(%) with positive coefficients will be 
introduced later. 



One-particle problem For j = 1, 2, ... and fixed a put 



.o,_ r o 17 .1.3...(2j-l)_ r o li (2j)! r o ■ 1 



ft = »(o) = U J -^f — - = ioYihtni ~ UV 



■2> L 2 J j\ L 2 J 2Jj\j\ L 2 J v^FJ 

Then we have 

Lemma 4 For j = 5, 6, ... the following inequalities hold 

> 1 [ S ] _ ^ fe+1 (£; + !)(*; + 2) 0jl o, fe 



J fc=l i 1+ ... + fe-m-l 2 2 Ji + 1 J fe + 1 

Proof. Let the particle, situated initially at the point x(0) = 0, move with the speed (a > 
being arbitrary) 

1 oo 

VI -at j=Q 

in the field of the external force F(x), which we should find. Then 

x(t) = / v(s)ds = ( — )Vl - at + -=>l-at = (l- ^-f 

Jq a a 2 

dv a 1 a 1 

F 



dt 2(1 -a*)! 2(l-f)3 
F<*> _ a, fc+1 3.4...(fc + 2) _ .a. fc+1 (fc + l)(fc + 2) 



fc! y 2 J k\ y 2 J 2 

Similarly to the above derivation of the recurrent equations (the only difference is that here 
v(0) = 1 and there are no A-terms), we get 

m = EE £ 7^ (fc) (^(o))^(o)-^-...-^- 

for Oj defined above. Here, as above, J2j 1 +...+j m =j- m -i i s the summation over all ji, ...,jk- P such 
that 

Ji + ••• +ifc- P =j-k-l 

As all coefficients are positive, then, neglecting the terms with ]? > 0, we have for any a > 

[— ] 

_ v^ y^ ( a \k+i 9ji 9j k ^ 

J fc=i ii+...+im=i-m-i 2 Ji + 1 Jfe + 1 

<Iy^ y- / a ^+i ( fc + 1 )( fe + 2 ) gn ft. < 

" J jfe Jl+ ... + ^- m -i 2 2 * + O* + 1 " 



Majorization From the recurrent formula for c^ one can see that they can be written as 

dij 
Cij = 2_^ Vi,j,r-L* l ' 3 ' r Mi,j,r 



where 6jj )7 . amd dij are some numbers not depending neither on N nor on F, and M it j^ r G A. 
By lemma [2] 

deg M i:jtr < j - 1 

We need other preliminary definitions. For any polynomial 

P = Y,b r N Ir M r 
where b r are the numbers, not depending neither on N nor on F, M r G A, we put 

Q(p\ = X^\b \]yIr+Q(M r ) £,Qo(M r )-Q(M r )+de S M r 
r 

where, for any monomial M r , the natural number Qo(M r ) equals the sum J2k over ai l factors 
(variables) F ijkjq . In particular, 

8(c-) = Y^ \b- ■ \]sr Ii -i- r+ Q( Mi >3>^c Qo ^ Mi ' J ' r ^ Q<KMi ' : '' r ' >+degMi ' j ' r 

r 

By definition for any i, k 

P{F itk>0 ) = C^ 1 ,/3(V ± F iifc , ) = p{F itk>1 ) = C k F +2 N~ l = C F N~ l p(F hk ) 
Moreover, for any two polynomials Pi, P 2 

P(Pi+P2) < /?(Pi) + /3(P 2 ),/?(PiP 2 ) < /?(Pi)/?(P 2 ) (24) 

and for any monomial M 

/3(V ± M) < (degM)N Q ^- 1 C? oiM) - QiM>HBgM+1 = (deg M^N^^M) (25) 

It follows that for any polynomial P 

f3(V ± P)<(degP)C F N- 1 f3(P) (26) 

We call f3(P) the majorant of P as by 

|*i,fc,i| = |V + F Jifc | <f' \F hk+1 (x)\dx < C k F +2 N~ l = 0(F iik>1 ) 

J Xi 

we have 

1^1 < P{P) 
From (El and (TT7D it follows that 



[^] [¥] 



0(cy)<E E P(Aij(rn;ji,...,j m ))+Y, E /3(B ij (k ]Jl ,...,j k )) 

m=\ ji+...+j m =j-m-l fc=l ji+-+jk=j-k-l 

Our inductive hypothesis will be (c gj = gj(l)) 

Pica) < X J N* j -§ gj J = 1,2,..., J - 2 (27) 



Initial data One can choose xo > so that for j — 1, 2, 3, 4 

xiN^-lgj > (3(c tJ ) 
In fact, only for j = 3 there is dependence on N, but |3 — | is exactly N. 

Inductive step for A-terms with m > 1 To estimate A-terms we distinguish two cases: 
m = 1 and m > 1. For m > 1 we use obvious bounds 



/3(V%) < 2/3(c u ),/?(V-(V + c iJl ...V + c JJm )) < 2 m+1 I]/?(Q 



j?/ 



p 



Then from flU} and ([20]) we get 



J 3\ + 1 Jm + 1 



p=l -?P 



<- ^ jy-2+mom+l J-m-1 ?y-|(J-m-l)-mj TT Vjy < 



< 2 m+l x J-m-l iv |j-f m + 1 ( ^ _g*E ' 

^ „=i 3p + 1 



as the exponent over iV for m > 2 has the bound 



5/t ^ 3 5 T 8 7 5 T 3 

2 + rriH — (J — m — 1) — m- = -J — m — I — < -J 

6 V ; 26 6 6 ~ 6 2 



Then by lemma H] with a = 2 (if % > 2) 



E E ^(^j(w;ii,...,j m )) < 

m = 2 ii+...+im=i— m— l 

3-1 

< iV^" V ± 2 m+i x— 1 e ^(fi-^rT)< 

o „-. i i -; : ™ 1 <^ 1 /ti ~T~ -L 



m=2 3i+... +3m,=j-m- 



-1 ^ p=l -?P + 



< (-) 3 JVt J -i x ^ (28) 

A 

Inductive step for A-terms with m = 1 In case m = 1 we shall use the following bounds 
for j = J - 2 > 3. By (J2SI), (I2ZD and lemma |2] 

/3(V + c,,) < (j - l))N- l C F p{ Clj ) < (j - l))N- 1 C FX J Nl^lg J 

Similarly 

/3(V-V + c M )| < ((j - l)C F ) 2 N-Wb-l 9j 

This gives additional summand (—2) in the exponent over N, equal to 



that is 



5/t n 3 5 T 3 

3-2 + - (3-2 - - < -J- - 

6 V ; 2 - 6 2 



A,(l; h) = Ml; J -2) = jld.lN ^ V j + -l~ 2 ~ 



< 2C FX >- 2 Nl J -l 9j _ 2 < 2 -^^X 1 Nl J -lg J (29) 



X 2 9 



Inductive step for 5-terms For S-terms the inductive bound is easier, but here the degree 
of monomials increases 

J K - J\ i" J- Jk "T 1 



1^ +1 l l R(c- ■ ) 0(c--)<- C * +1 9jl -^ 



< - F . . . -— ^(c,;,A../3(c^O < ~ F . .^ : ...-Mh-yJi+—+iu N h+—+J 



k 



By lemma H] 



I 2 ^] J IS 1 ] £*+l 









3 fe=i ii+...+ik=i-fc-i -?' 1 + - 1 ^' fc + x 



(30) 



Finally, to end the proof of the theorem, we sum up three obtained expression (128|) . (129]) . ( 130]) . 
and choose x = Xi > so that 

(Sx- 1 + ^^ + C F e^) X - 2 < 1 

Then for any x > max(xo) Xi) we wn l have 

\cij\ < X J Ni gj < X J Ni 
The theorem is proved. 
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